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abstract 

We consider compactification of extra dimensions and numerically calculate Casimir 
energy which is provided by the mass of Kaluza-Klein modes. For the extra space we 
consider a torus with shape moduli and show that the corresponding vacuum energy is 
represented as a function of the moduli parameter of the extra dimensions. By assuming 
that the Casimir energy may be identified with cosmological constant, we evaluate the 
size of extra dimensions in terms of the recent data given by the Wilkinson Microwave 
Anisotropy Probe (WMAP) measurement and the supernovae observations. We suggest 
that the observed cosmological constant may probe the shape moduli of the extra space 
by the study of the Casimir energy of the compactified extra dimensions. 
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1. Introduction 

Extra dimensions have been introduced in order to solve the hierarchy problem |]J, 
and a lot of related works including the milestone papers and later developments 

have been done. In these approaches, gauge fields are localized on a four dimensional 
brane, which is regarded as our real world, and only a graviton can propagate in the extra 
space, which is transverse to the brane. When the extra space is compact, there exist 
the Kaluza-Klein modes of the graviton. They contribute to Casimir energy on the four 
dimensional brane. 

On the other hand, the origin of cosmological constant is also a standing problem. 



Recently the supernovae [|IO| , |TlH and WMAP |L2l , |l3H observations have revealed some pa- 



rameters of our universe. These new results encourage not only the researchers of cos- 



mology but also the particle physicists [14-16| toward challenging the resolution of the 



unsolved problems again. If we assume that the Casimir energy, or in other words the 
vacuum energy, can be regarded as the cosmological constant, we can estimate the size of 



extra dimensions in terms of the observational value of the cosmological constant Jl6j . 

When we calculate the vacuum energy, we have to perform a summation over the 
infinite number of Kaluza-Klein modes. In this paper, we introduce a new mathematical 
technique to achieve this infinite summation, which provides a reasonable regularization 
and makes it possible to calculate numerically the vacuum energy of all Kaluza-Klein 
modes. 

In it is shown that the number of extra dimensions should be two from the 

consistency with the hierarchy problem. The choice of the two extra dimensions are also 



favourable from the viewpoint of the cosmological constant problem [IE|. So we consider 
the extra space of two dimensions in this paper. In particular we concentrate on a torus 
with shape moduli. The compactification on the torus and the Kaluza-Klein masses under 
a standard mass operator are analysed in |fTT|-p2||. In this paper we additionally impose 
(Z2) 2 symmetry or Dirichlet conditions on the torus. Since the standard mass operator is 
inconsistent with these boundary conditions, we propose a new mass operator incorporating 
the shape moduli in consistency with the boundary conditions on the torus. Using this 
operator, we calculate the vacuum energy. The resulting expression depends on the shape 
moduli parameter, as a function of which the structure of the extra space can be probed 
by confronting the numerical values of the Casimir energy with the observed cosmological 



constant 10-13 . 



1 



The paper is organized as follows. In the section 2, we consider the (Z 2 ) 2 symmetry 
and Dirichlet conditions on the torus with the shape moduli. Requiring the invariance 
under these conditions, we introduce the new mass operator with the shape moduli. In 
the section 3, we show the technical details of the regularization of the vacuum energy 
to which the infinite number of Kaluza-Klein modes contribute. We also present the 
numerical results of the calculation of the vacuum energy. In the section 4, we identify 
these numerical values of the vacuum energy with the cosmological constant which is given 
by [pT0| pT3|1 and evaluate the size of extra dimensions. The section 5 is devoted to conclusions 
and discussion. 



2. Shape Moduli and Mass Operator 



Though various kinds of extra space are available, we consider a two dimensional 
torus, T 2 , in this paper. In the covering space of T 2 , we require in the standard way the 
two periodic conditions, 



(2/1,2/2) ~ (2/1 + 27r-Ri,y 2 ), 

(2/1,2/2) ~ (2/1 + 2nR 2 cos 9, y 2 + 2irR 2 sm9). 



(2.1) 



Ri and R 2 denote the radii of two cycles, while 9 parametrizes the shape moduli of T 2 . 
After the Kaluze-Klein compactification on such a torus, the eigen functions for the usual 
(mass) 2 operator — [(d/dyi) 2 + (d/dy 2 ) 2 ] are proportional to 



exp 



. ni 

'777 1 



2/2 

tan 9 



+ 1 



. n 2 y 2 



R 2 sin 9 



(2.2) 



where ni,n 2 G Z, and the eigen values become n\/ R\+n 2 2 / R 2 2 -2n \n 2 cos9/(RiR 2 ). 

For these mass eigen values, we can numerically calculate, in the same way as in ||16|1 , the 
vacuum energy in which we have to sum up all the Kaluza-Klein modes. 
Now we introduce new coordinates (Xi,X 2 ) defined by 



X X EE 



2/1 



2/2 



Ri l?itan# 
The periodic condition ( |2.1|) is rewritten as 



Xo EE 



2/2 



R 2 sin 9 



(2.3) 



{X 1 ,X 2 )~(X 1 + 2tt,X 2 ), 
(X 1 ,X 2 )~(X 1 ,X 2 + 2n). 



(2-4) 
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The eigen functions (|2.2|) are equal to exp(mi-Xi + ^2X2) with ^ £ Z (i = 1,2). 
In addition, let us impose (Z2) 2 symmetry on T 2 . 




>yi 



2%Rr 

fig. 1: A fundamental region of T 2 in the covering space. 

This symmetry is realized by the two mirror symmetry along the two thick lines represented 
in fig. 1. In terms of the coordinates (2/1,2/2)5 the symmetry is denoted by the following 
identification 

(2/1,2/2) 

\ 1 mm n 1 

(2 - 5) 

(2/1,2/2) ~ ( yi - 7—^77, -2/2 



22/2 

"fl + I a^ 2 

tanc/ 

22/2 



tan#' 



In fig. 1, the point A'(— y\ + is identified with the point A(yi,y 2 ) by one Z 2 

symmetry and the point A" (2/1 — tl^, — 2/2) is also identified with the point A by the 

^2 



other Z2 symmetry. Since the (Z2) symmetry ( j2.5|) can be interpreted as the reflection 
symmetry given by 

(X 1 ,X 2 ) ~ (-A l7 A 2 ), (A l5 A 2 ) ~ (A l7 -A 2 ), (2.6) 

the wave function under the (Z2) 2 symmetry, through the use of the coordinates (Xi, X2), 
should be an even function with respect to each variable, X\ or X%- By the analogy of ( |2.2|) , 
we may suppose that the wave function satisfying ( |2.1|) and ( p.5|) is cos(niAi) cos(n2A2). 
But this function is not an eigen function of the (mass) 2 operator — [(d/dyi) 2 + (d/dy2) 2 ] ■ 
So we propose a new mass operator which has cos(ni-Xi) cos(n2A2) as an eigen func- 
tion. First we note that the following differential operators simply change signs under the 
operation of (Z2) 2 symmetry on T 2 , ( |2.5| ): 
did d 



d yi R x dX x 



+ 



1 d 



1 



d 



dy2 tan 6 dyi R2 sin 6 8X 2 



(2.7) 
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Though there are a variety of choices for the mass operator, the above observation has led 
us to consider, as the simplest toy model, the following generic (mass) 2 operator defined 
by 



1 d_ 



+ 



1 



d 



R 2 sin9dX 2 



1 d 2 d 2 2 d d 

+ 7^7 + 



sin 6 dyf dy\ tan 9 dyi dy 2 



[2.1 



The proposed (mass) 2 operator ( p78|) is invariant under the (Z2) 2 symmetry (|2.5|), and 
reduces to the standard one for 9 = 7r/2, that is, when the fundamental region of T 2 is 
rectangular. 

In other words, though we can make a variety of choices for the metric in the extra 
space, the (Z 2 ) 2 symmetry on T 2 has guided us to choose the simplest case of (|2.8|) for the 
mass operator. The corresponding metric of T 2 can be proved to be 



dyi- 



tan 9 



(dy 1 dy 2 + dy 2 dy x ) + 



sin 9 



This metric is naturally invariant under the (Z2) 2 symmetry ( |2.5| ), and again reduces to 
the standard rectangular one when 9 = 7r/2. 

In the next section, we use the mass operator ( |2.8| ) and calculate the vacuum energy. 



3. Vacuum Energy and Boundary Conditions 

We consider the torus compactification with shape moduli. The torus T 2 without 
shape moduli is the same as {S 1 )® 2 . The vacuum (Casimir) energy which is enhanced by 
the Kaluza-Klein modes derived from the compactification on the d dimensional extra space 
(S l )® d has already been discussed in | l6[j . In this section we consider the compactification 



on T 2 with shape moduli and the generic mass operator (|2.8| ) as has been discussed in the 
section [2| We then calculate the vacuum energy in our four dimensional space to which 
all the Kaluza-Klein modes from the extra space contribute. 
For simplicity we set Ri = R 2 = R. 

3.1. without a boundary condition 

The eigen functions compactified on T 2 with the mass operator (|2.8j) are described by 



exp(miXi) exp(in2^2) from (|2.4j), and the masses of Kaluza-Klein modes then become 



1 / n 2 1 

M = ~R \ n ? + "Hb = p/m,n a ((9), n u n 2 e Z. 
K V sin 9 ri 



We consider the zero point energies of these Kaluze-Klein modes and the vacuum 
energy of our four dimensional world is represented by 



E 



1 oo 

\ E 



ni ,n2 = — oo 



d 3 k 

(2tt) 3 



k 2 + M 2 . 



(3.1) 



Note that, if the Kaluza-Klein modes are for gravitons, we have to insert a factor cor- 
responding to the number of polarization degrees of freedom, (2 + 2)(3 + 2)/2 — 1 = 9, 
in front of the right hand side of ( |3.1| ) | |23| |. In the following we consider only a scalar 
degree of freedom for simplicity. Since the integration in ( |3.1| ) is generally divergent, we 
use dimensional regularization for it and obtain 



E 



1 



oo 



64tt 2 ^ 

TL\ ,712 = — OO 



M z 



(3.2) 



where the parameter A has a mass dimension. Since the infinite summation for 
n\ + n|/ sin 2 9 J is calculated as 



oo oo s 1 \ oo 

E W) 4 = E ( 1 + ^K+ E 



Til i"-2 = — OO 



n,m= — oo 



n,m= — oo 



sin 6 9 



2 2 

-n m 



= 4(1 + 
= 0. 



C(0)C(-4) + t^C(-2) 



sin 4 # / sin 2 # 



where ((— 2) = £( — 4) = 0, the vacuum energy ( |3.2|) is reduced to 



E 



1 



32tt 2 j R 4 
1 



77 2 



Til ,77,2= — OO 
OO 



32tt 2 j R 4 ^ dx' 

Til ,77,2 = — OO 



11 2 



x=4 



(3.3) 



In calculating the vacuum energy, the infinite summation shows up and appears to diverge 
because of the contribution of the infinite number of Kaluza-Klein modes. It is known 
that in some special cases [|4[^5[ the vacuum energy can be regularized and are calcu- 
lated exactly in terms of Epstein's zeta function. Since in general we can not use such a 



regularization, we suggest a new numerical method as follows. By the use of the Fourier 



transformation [26 



2-sin4f / x\ „ fx + d\ f°° expfeti^ 



7T2 



(3.4) 



which is valid for any real values of (a±, 0*2, • • • , cr^), we obtain 

00 00 

fn 1 ,n 2 = 



2 a; sin^ r / | x 



7T 



2 r 00 exp(imgi+i^gf 2 ) 



m ,n2=—oo 



n\ ,n2 =—oo 



x+2 



sin 



7T 



2+x 



771 1 > m 2 = — 00 
( mi ,m 2 )^(0,0) 



(mf + m,2 sin 0) 



(3.5) 



where we have used ^ n= _ 00 exp(m<z) = 2n ^ m= _ oo 5(5 — 2itm). The contribution of 
(mi, 7712) = (0, 0) mode vanishes, which can be proven by employing a generalized function 
technique |2~6|| . From (|3.3| ) and (|3.5|), the vacuum energy becomes 



E 



64n 7 R 4 



^ 00 
T - i r(3) 2 sin0 £ 



m 1 > m 2 = — 00 
(m 1 ,m 2 )^(0,0) 



(mf + m 2 sin 2 6*) 3 



(3.6) 



Since ( |3.6|) is convergent, next we numerically calculate it. 
R 4 E 



0.25/"0.5 0.75 1 1.25 1.5 



6 [rad] 



-0 . 02 

-0 . 04 
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-0.08 
-0.1 
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fig. 2: The 9 dependence of the vacuum energy with no boundary condition. 

The values are negative for all 9. 



The numerical value of R 4 E which depends on 9 is represented in fig. 2. R 4 E is negative 
definite for all 9. 
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3.2. with (Z2) 2 symmetry 

We consider (Z2) 2 symmetry on T 2 . Since the wave function should be an even 
function of X% and X 2 respectively on account of the condition (|2.6|), the eigen function 
for the mass operator (|2.8[) becomes proportional to 



cos(niAij cos(n2A 2 J = costt-i ( — ^— ] cosn2 



We then obtain the masses of Kaluza-Klein modes, 



R2 sin 9 



m,n 2 e {0,N}. 



M=-f ni ,n 2 (0), m,n 2 6{0,N}. 



(3.7) 



Note that, since the eigen functions consist of cosine functions, ni and n 2 are zero or 
positive integers. In terms of ( j377l) , we analyse the vacuum energy, and obtain 



64tt 



1 00 

h £ ^ 



ni,ri2=0 
oc 



\ A 2 3 
In — ^ + - 

M 2 2 



.2 p4 ^ ] fn 1 ,n 2 hl/ ni 



327T 2 i? 4 

1 

32^R4 
1 

128^R4 



"2 



ni,«2=0 
00 



III ,Tl2 = — O0 
OO 



n2 = — 00 



,"•2 



- ^ (1 + ^,0) 2 E ( 1 + ( W) 

ni = — 00 J L 

00 

i 2 =— 00 
00 00 



/n 1,712 l n /tii 



"2 



"2 



n\ = — 00 



ri2=— 00 



(3. 



where the last term has been dropped due to the identity /o,o 4 hi/0,0 = 0. The first term 
in ( ^.8| ) has already been considered in the previous subsection, and from (|3.5| ) we get 

d 



} ^ /ni,H2 l n /ni,n 2 ^ ^ dx^™ 1 '™ 2 



ni ,n 2 =— 00 



ni ,n 2 = — 00 



x=4 
00 



L r( 3) 2 sin0 J2 



2tt 5 



771 1 i m 2~ — 00 
(m 1 ,m 2 )#(0,Q) 



(m 2 + m 2 , sin 2 #)' 



(3.9) 



Using (|3.4|) , we calculate the second term in (|3.8|) as 



Til = — OO 



ni = — oo 



x=4 



3 _l_4 



V 

^ r(3)r V 2 

27T2 V 1 



-x-1 



\ 



k = — oo 



OO 

£ IT 

fc^O 



4-r(3)r(|)C(5). 



x=4 



(3.10) 



We can also show in the same way as for (|3.10| ) that the third term in (|3~8|) becomes 



/o,n 2 4 ln/ ,n 2 --^ Yl 



a-, 



ri2 = — oo 



ri2=— oo 



sin 2 # 



x=4 



dx \ (sm9) x ^ 



ri2=— oo 



x=4 



1 



, „-r(3)r ( - ) c(5). 



7T2 sin 6 



(3.11) 



Substituting (|3.9|), ( |3.10| ) and ( |3.11| ) into ( |3.8|) , the vacuum energy with the condition of 
(Z2) 2 symmetry is written down as 



%2 256yr 7 i? 4 



k,l = — 00 

(M)#(o,o) 



r(3) 



2567r 7 i? 4 



4r(3)sin# £ 



+ 2T(3)C(8) (sin, + + 2^ r g) C (5) (l + 



(3.12) 



8 
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fig. 3: The 9 dependence of the vacuum energy with (Z2) 2 symmetry. 
The values are negative for all 9. 

We can calculate the coefficient of R~ 4 in ( |3.12| ) numerically and the result is depicted in 
fig. 3. The vacuum energy E% 2 is negative for all 9. 

3.3. with Dirichlet condition 

We require the Dirichlet condition on the boundary of the fundamental region of T 2 , 
that is, on the two lines defined by y± tan 9 — y% = and 7/2 = 0. In other words, the wave 
function ^(X±, X 2 ) should vanish on the lines X\ = and X 2 = 0. So the eigen functions 
of ( |2.8| ) are proportional to the products of sine functions, 



sin(niXi) sin(n2X2) = sinni 



Ui 



y-2 



SUl 77-2 



U2 



Ri i?itan6* / \R2sm.9 
The eigen values, that is, the masses of Kaluza-Klein modes become 

1 



ni, 7T-2 G N. 



M 



R 



f ni ,n 2 (9), ni.njGN. 



(3.13) 



By the use of ( |3.9[ ) , ( |3.1U[ ) and ( |3.11| ) , we obtain the vacuum energy, 



-Ebiri 



chlet 



64tt 2 
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32n 2 R 4 



b £ * 



ni,n,2 = l 
00 



In — 7T + - 

M 2 2 



/rei,ri2 l n /ni, 



??, 2 



ni = — 00 



n2 = — 00 



9 



1 



128tt 2 j R 4 

r(3) 



256n 7 R 4 



ni = — QG 
1 



n.2=— oo 



^ /ni,n 2 4 l n /ni,n 2 — /«i ,0 4 m /m ,0 ~ /o,n 2 4 l n /o 

,ni ,7i2 = — oo 

oo 

4r(3)sin# V 

/K\ / 1 \ 

(3.14) 



^ (n? + n 2 , sin 2 9) 3 

ii,ri2=l v 1 z ' 

+ 2r (3 )c(6) (sin* + -^) - ^ r (§) C (5) (i + -L_) 



R £ n : 



Dirichlet 
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(b) 





Tl e[r 



fig. 4: The 9 dependence of the vacuum energy with Dirichlet condition. The values turn 
from negative to positive at 9 = 9 CT as shown in the magnified plot (b). 

The results of the numerical calculations of -R 4 -E , Dirichiet are represented in fig. 4. In 
this case it is significant that the positive vacuum energy appears. The vacuum energy 
^Dirichlet is equal to zero at 9 = 9 CT (fig. 4(b)). For 9 < 9 cr the vacuum energy is negative, 
while for 9 > 9 CT it is positive. The critical angle 9 cr is approximately equal to 0.947243 
[rad] which is calculated numerically. 

The positive vacuum energy is favourable, because the cosmological constant is posi- 
tive. From this point of view, we will give some comments in the next section. 



4. Comments on Generic Size of Extra Dimensions 

The extra dimensions have been introduced in Jl]||] in order to solve the hierarchy 
problem . When we consider a 4 + d dimensional space-time, that is, a d dimensional extra 
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space, the relation between the Plank scale M p \ of our four dimensional world and the one 
Mpi( 4 _|_d) of the 4 + d dimensional space-time is described by 



M pl 2 ~ M m+d) 2+d R d , 

where R denotes the size of the extra dimensions. When we assume that M v \u +( j\ ~ 
toew ~ 1 [TeV] , we obtain the size for the extra space as 

„ 30 17r , 

R ~ 10~" 17 [cm]. 

For d — 2, the size R becomes a sub-millimeter size, 10 _2 [cm]. 

On the other hand, we have not had a complete answer for the cosmological constant 
problem yet. But a lot of features have been made clear by some observations, for a 
recent example, the WMAP observation. The results of the WMAP HT2|JT^1 provide us the 
following data; 

h = 0.71 (Hubble constant), (4.1a) 
fitot = 1-02 ± 0.02 (total density), (4.16) 
fi A = 0.73 ± 0.04 (dark energy density). (4.1c) 



From ( ^.la[ ), the Hubble constant is 

H = hx 100[km • sec -1 ■ Mpc" 1 ] = 0.7675 x 10" 28 [cm _1 ]. (4.2) 



In terms of ( fL2|) , we obtain the critical density, 

o tt2 

p c = — ^ = 5.313 x 10 3 [eV- cm" 3 ]. (4.3) 
87rG 

Since we consider that the dark energy is almost the same as the cosmological constant A, 
it is calculated from (|4.16p, (|4.1c[) and (|4.3|) as 



A = p c O to tfiA = 3.956 x 10 3 [eV ■ cm -3 ]. (4.4) 

Now we assume that the cosmological constant is identified with the vacuum energy to 
which the Kaluza-Klein modes from the compact ificat ion of extra dimensions contribute. 
The vacuum energy is proportional to R . If the size of the extra space R is a submil- 
limeter size, i? -4 is in the order of 10 8 [cm -4 ] ~ 10 4 [eV • cm -3 ]. Since it is close to the 
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cosmological constant (|4.4|) , that identification is reasonable. The relationship between 
the vacuum energy, in other words, the Casimir energy, and the cosmological constant has 



been discussed in the various papers []T6|J2~7| - |37|1 . We believe from recent observations that 
the cosmological constant is positive. On account of that identification, the vacuum energy 
also has to be positive. Fortunately we were able to realize a positive vacuum energy by 



the use of the extra dimensions with the Dirichlet condition in the section 3.3. 



^Dirichlet 




fig. 5: The R dependence of the vacuum energy. If \9 — 9 cr \ increases, the corresponding line 
for the graph of -^Dirichlet) moves away from the vertical axis. 



Since the vacuum energy -Ebmchiet depends on R~ 4 from (|3.14|) , the behaviour of .Ebirichiet 
for R is represented by fig. 5. 

If the values of the cosmological constant A and the moduli parameter 9 are given, the 
size of extra dimensions R is determined. For example, let us consider the case of 9 = tt/2. 
This is the situation in which the torus T 2 has no more shape moduli and is essentially equal 
to the product of two spheres, (5 1 )® 2 , with the same radius R. The vacuum energy ( |3.14 ) 
then becomes £ D irichiet = 1-1769 x l(T 6 /# 4 - In terms of he = 1.973269 x 10" 5 [eV • cm], 
we obtain 

2.3223 x 10" 11 

^Dirichlet eV - cm = — — r— . 

(it[cmj) 4 

When we identify this vacuum energy with the cosmological constant Q4.4| ) , the size of the 
extra dimensions R becomes 2.768 x 10 _4 [cm]. Though this scale is different by 10 -2 order 
from the one derived from the hierarchy problem [l],[| , it is not so crucial an inconsistency. 

Note that there is no zero mode under the Dirichlet condition. This is a serious 
problem, if the particle propagating in the extra space is a graviton. But it is not a 
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problem for a scalar particle, because, in fact, there is no such a massless scalar particle 
in our four dimensional space-time. 



5. Discussions and Conclusions 

We considered the torus T 2 as a two dimensional extra space. We introduced the shape 
moduli of T 2 and analysed the dependence of the vacuum energy, or the Casimir energy, 
on the shape moduli parameter, to which all the Kaluza-Klein modes derived from the T 2 
compactification contributed. Without a condition on the covering space of T 2 except for 
bi-periodic condition, we can find eigen functions for the standard mass operator. But, if 
we require the (Z2) 2 or the Dirichlet conditions, there is no eigen function for the mass 
operator. So we proposed the new mass operator (|2.8j ) together with the corresponding 
eigen functions and considered the vacuum energy for this operator. 

In terms of the mathematical tools of the Fourier transformation and the generalized 
functions, we regularized the divergence of the infinite summations of Kaluza-Klein modes. 
For a specific case of d = 1, our regularization reduces to the zeta- function regularization, 
where higher mode is suppressed and infinity arising in infinite volume limit is subtracted. 



This regularization provides us the same results as are given by the cuttoff method ||38|| . 
Our method being generalized to multi-dimensional problems is a valid and unique one to 
the extent that the zeta-function regularization is a valid method. 

We then made it possible to calculate numerically the vacuum energy including the 
contributions of all the Kaluza-Klein modes. In general the vacuum energy has a negative 
value. But we have found that, if we impose the Dirichlet conditions, the value of the 
vacuum energy turns out to be positive lifting from a negative territory, as the shape 
moduli parameter changes crossing the critical value. 

We assumed that the vacuum energy was identified with the cosmological constant. 
Since it is confirmed by the observations, for example, the WMAP observation, that the 
cosmological constant is positive, we compared the WMAP's observational value of the 
cosmological constant with the vacuum energy which we calculated under the Dirichlet 
condition. From this comparison we showed that the size of extra dimensions is of order 
10 _4 [cm]. The value of R = 2.768 x 10 _4 [cm] presented in the section 4 is the maximal 
one that we can realize in the model considered in this paper. On the other hand, from the 
hierarchy problem, the size of extra space for the two extra dimensions is 10 ~ 2 [cm] . Though 
these two sizes do not accurately agree with each other, they are sufficiently consistent. 
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Since the cosmological constant may consist not only of the vacuum energy but also 
of some other elements, it is significant that we showed that the two sizes of extra dimen- 
sions derived from the hierarchy problem and from the cosmological constant problem are 
similar. 

It is possible to consider various models with an extra space. So one of future problems 
is for us to find the model in which the vacuum energy has a positive value of right order 
of magnitude. Since we calculated the vacuum energy only in a semi-classical limit in this 
paper, it is also a future problem to consider possible quantum loop effects and also some 
possible contributions coming from quantum fluctuations of relevant metric components. 
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